State-of-the-art stellarator optimization code
This tutorial walks the user through running TERPSICHORE on an NCSX equilibrium.
TERPSICHORE in standalone mode requires the user to provide a text file (named fort.18) which described the VMEC equilibrium. One can run STELLOPT (with kink targeting turned on) to generate these files (terpsichore_eq.ext files). In this example we provide one based on the LI383 equilibrium. examples/fort.18
The TERPSICHORE source must be modified for a given set of parameters. The file to edit is tpr_modules_ap.f. Edit the header as so
NI = 48
NJ=151
NK=151
MLMNV=94
MLMNB=256
MLMNS=382
IVAC=24
LSSL=4446
LSSD=2745
MMAXDF=127
NMAXDF=31
The simply rebuild the code by running the make command.
With the code compiled you'll need to generate an input file. The one below can be used to examine the n=1 family of modes in the equilibrium. You'll notice in the stability mode table that not just the n=1 modes are examined. This is because in stellarator geometry the modes couple by +/-NFP so for NFP=3 we have n=...-5,-2,1,4,7...
QAS-3n0
C
C MM NMIN NMAX MMS NSMIN NSMAX NPROCS INSOL
17 -6 +8 25 -20 22 1 0
C
C TABLE OF FOURIER COEFFIENTS FOR BOOZER COORDINATES
C EQUILIBRIUM SETTINGS ARE COMPUTED FROM FIT/VMEC
C
C M= 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 N
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -6
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -5
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -4
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -3
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -2
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 5
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 6
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 7
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 8
C
LLAMPR LVMTPR LMETPR LFOUPR
0 0 0 0
LLHSPR LRHSPR LEIGPR LEFCPR
9 9 1 1
LXYZPR LIOTPL LDW2PL LEFCPL
0 1 1 1
LCURRF LMESHP LMESHV LITERS
1 1 2 1
LXYZPL LEFPLS LEQVPL LPRESS
1 1 0 0
C
C PVAC PARFAC QONAX QN DSVAC QVAC NOWALL
1.0001e+00 0.0000e-00 0.6500e-00 0.0000e-00 1.0000e-00 1.0001e+00 -2
C
C AWALL EWALL DWALL GWALL DRWAL DZWAL NPWALL
1.2000e+00 2.0000e+00 -0.3200e-00 5.2000e-00 0.0000e-00 -0.1000e-00 3
C
C RPLMIN XPLO DELTAJP WCT CURFAC
1.0000e-05 1.0000e-06 4.0000e-02 6.6667e-01 1.0000e-00
C
C MODELK = 0
C
C NUMBER OF EQUILIBRIUM FIELD PERIODS PER STABILITY PERIOD: NSTA = 3
C
C TABLE OF FOURIER COEFFIENTS FOR STABILITY DISPLACEMENTS
C
C M= 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 N
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1-20
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-19
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-18
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1-17
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-16
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-15
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1-14
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-13
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-12
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1-11
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-10
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -9
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -8
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -7
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -5
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -4
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -3
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -2
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 6
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 7
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 8
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 9
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 10
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 11
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 13
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 14
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 15
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 16
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 17
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 18
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 19
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 21
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 22
C
C NEV NITMAX AL0 EPSPAM IGREEN MPINIT
1 4500 -0.500E-04 1.000E-04 0 0
C
Now we can execute the code by passing the TERPSICHORE input file to the code through standard in
>> ./tpr_ap.x < terpsichore_input_01
First some information from the input files is printed on the screen.
1
###### TERPSICHORE ######### EQINVM ######### TERPSICHORE ####
QAS-3n0
NIM IVAC NJ NK MM NMIN NMAX
48 12 421 423 17 -6 8
TABLE OF R/Z (PR/Z) COEFFICIENTS
M= 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 ...N
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -6
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -5
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -4
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -3
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -2
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 3
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 5
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 6
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 7
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 8
PRESSURE GRADIENT RECONSTRUCTED FROM RADIAL FORCE BALANCE
PARALLEL CURRENT COMPUTED FROM CHARGE CONSERVATION
QN QONAX AIOTAX PARFAC SVMAX - 1 PVAC QVAC NOWALL
0.0 0.6500 1.5385 0.0000 1.000E+00 1.000E+00 1.000E+00 -2
npwall awall ewall dwall gwall drwal dzwal
3 1.200E+00 2.000E+00-3.200E-01 5.200E+00 0.000E+00-1.000E-01
rplmin = 1.0000E-05 xplo = 1.0000E-06 djp = 4.0000E-02 wct = 6.6667E-01 curfac= 1.0000E+00
The equilibrium profiles are output VMEC (iota, plasma mass, differential volume, pressure)
IOTA M VP P
0.00E+00 0.00E+00 0.00000000E+00 0.00E+00
2.25E-01 4.47E+04 2.25268473E+00 5.61E-02
2.37E-01 4.46E+04 2.24866205E+00 5.61E-02
2.48E-01 4.45E+04 2.24454039E+00 5.60E-02
2.58E-01 4.44E+04 2.24029036E+00 5.58E-02
2.67E-01 4.42E+04 2.23588561E+00 5.56E-02
2.75E-01 4.40E+04 2.23129364E+00 5.53E-02
2.83E-01 4.37E+04 2.22648029E+00 5.49E-02
2.91E-01 4.33E+04 2.22141336E+00 5.44E-02
2.99E-01 4.29E+04 2.21606545E+00 5.39E-02
3.06E-01 4.24E+04 2.21041636E+00 5.33E-02
3.13E-01 4.19E+04 2.20445449E+00 5.26E-02
3.19E-01 4.13E+04 2.19817763E+00 5.18E-02
3.26E-01 4.06E+04 2.19159288E+00 5.10E-02
3.32E-01 3.99E+04 2.18471611E+00 5.01E-02
3.38E-01 3.91E+04 2.17757078E+00 4.91E-02
3.44E-01 3.82E+04 2.17018675E+00 4.80E-02
3.50E-01 3.73E+04 2.16259857E+00 4.69E-02
3.56E-01 3.63E+04 2.15484411E+00 4.56E-02
3.62E-01 3.53E+04 2.14696298E+00 4.43E-02
3.68E-01 3.42E+04 2.13899526E+00 4.30E-02
3.74E-01 3.31E+04 2.13098038E+00 4.15E-02
3.81E-01 3.19E+04 2.12295623E+00 4.01E-02
3.87E-01 3.06E+04 2.11495837E+00 3.85E-02
3.94E-01 2.94E+04 2.10701960E+00 3.69E-02
4.00E-01 2.80E+04 2.09916945E+00 3.52E-02
4.07E-01 2.67E+04 2.09143396E+00 3.35E-02
4.14E-01 2.53E+04 2.08383547E+00 3.18E-02
4.21E-01 2.39E+04 2.07639249E+00 3.00E-02
4.29E-01 2.24E+04 2.06911975E+00 2.82E-02
4.37E-01 2.10E+04 2.06202819E+00 2.64E-02
4.44E-01 1.95E+04 2.05512516E+00 2.45E-02
4.52E-01 1.80E+04 2.04841459E+00 2.27E-02
4.59E-01 1.65E+04 2.04189723E+00 2.08E-02
4.67E-01 1.51E+04 2.03557088E+00 1.89E-02
4.74E-01 1.36E+04 2.02943075E+00 1.71E-02
4.81E-01 1.21E+04 2.02346973E+00 1.53E-02
4.88E-01 1.07E+04 2.01767893E+00 1.35E-02
4.94E-01 9.33E+03 2.01204812E+00 1.17E-02
4.99E-01 7.99E+03 2.00656660E+00 1.00E-02
5.04E-01 6.70E+03 2.00122409E+00 8.42E-03
5.08E-01 5.48E+03 1.99601234E+00 6.89E-03
5.10E-01 4.34E+03 1.99092712E+00 5.46E-03
5.12E-01 3.29E+03 1.98597135E+00 4.14E-03
5.12E-01 2.35E+03 1.98115952E+00 2.95E-03
5.11E-01 1.53E+03 1.97652510E+00 1.92E-03
5.08E-01 8.54E+02 1.97212915E+00 1.07E-03
5.05E-01 3.49E+02 1.96809273E+00 4.38E-04
5.01E-01 4.94E+01 1.96456631E+00 6.20E-05
EPS 1.920E-01 EL 1.091E+01 NPER 3 GAM 0.000E+00 ALPHA 0.00E+00
LMNV 94 LMNL263 LMNB264
Number of Modes in the Lambda Reconstruction Table:1537
The recalculated equilibrium profiles are now output. The differential volume (VVP) should be the same as before. The residual force balance (EQUI) should also be less than 1E-2, otherwise the boozer modes should be increased (except perhaps at the edge)
***** LGIKVM, RECONSTRUCTED EQUILIBRIUM *****
I VVP(I) PVP(I) PVPI(I) PARPVI(I) CJVP(I) CIVP(I) EQUI
1 2.2526847E+00 4.1877E-03-2.1561E-03-0.0000E+00-7.8124E-03-2.0656E-02 8.8773E-01
2 2.2486621E+00-4.3983E-03-4.9130E-03-0.0000E+00-1.2256E-02 1.6792E-02 5.5269E-02
3 2.2445404E+00-7.7679E-03-7.9723E-03-0.0000E+00-1.5910E-02 3.0818E-02 1.2988E-02
4 2.2402904E+00-1.1133E-02-1.1236E-02-0.0000E+00-1.9107E-02 4.4823E-02 4.6174E-03
5 2.2358856E+00-1.4594E-02-1.4651E-02-0.0000E+00-2.1963E-02 5.9245E-02 1.9337E-03
6 2.2312936E+00-1.8151E-02-1.8182E-02-0.0000E+00-2.4708E-02 7.4014E-02 8.3407E-04
7 2.2264803E+00-2.1787E-02-2.1801E-02-0.0000E+00-2.7320E-02 8.9059E-02 3.2322E-04
8 2.2214134E+00-2.5482E-02-2.5485E-02-0.0000E+00-2.9938E-02 1.0425E-01 6.6410E-05
9 2.2160654E+00-2.9220E-02-2.9216E-02-0.0000E+00-3.2487E-02 1.1953E-01-6.1021E-05
10 2.2104164E+00-3.2982E-02-3.2974E-02-0.0000E+00-3.5083E-02 1.3476E-01-1.1930E-04
11 2.2044545E+00-3.6753E-02-3.6743E-02-0.0000E+00-3.7641E-02 1.4992E-01-1.4212E-04
12 2.1981776E+00-4.0516E-02-4.0505E-02-0.0000E+00-4.0277E-02 1.6486E-01-1.3646E-04
13 2.1915929E+00-4.4257E-02-4.4246E-02-0.0000E+00-4.2913E-02 1.7956E-01-1.1467E-04
14 2.1847161E+00-4.7959E-02-4.7950E-02-0.0000E+00-4.5667E-02 1.9390E-01-9.1167E-05
15 2.1775708E+00-5.1608E-02-5.1602E-02-0.0000E+00-4.8468E-02 2.0783E-01-6.3163E-05
16 2.1701867E+00-5.5190E-02-5.5185E-02-0.0000E+00-5.1432E-02 2.2125E-01-4.1436E-05
17 2.1625986E+00-5.8687E-02-5.8685E-02-0.0000E+00-5.4484E-02 2.3412E-01-2.2653E-05
18 2.1548441E+00-6.2086E-02-6.2085E-02-0.0000E+00-5.7732E-02 2.4633E-01-5.9264E-06
19 2.1469630E+00-6.5369E-02-6.5370E-02-0.0000E+00-6.1093E-02 2.5785E-01 5.9667E-06
20 2.1389953E+00-6.8522E-02-6.8523E-02-0.0000E+00-6.4659E-02 2.6857E-01 1.2819E-05
21 2.1309804E+00-7.1526E-02-7.1529E-02-0.0000E+00-6.8334E-02 2.7847E-01 1.9523E-05
22 2.1229562E+00-7.4366E-02-7.4370E-02-0.0000E+00-7.2191E-02 2.8745E-01 2.2693E-05
23 2.1149584E+00-7.7024E-02-7.7028E-02-0.0000E+00-7.6122E-02 2.9549E-01 2.4507E-05
24 2.1070196E+00-7.9482E-02-7.9487E-02-0.0000E+00-8.0180E-02 3.0250E-01 2.6896E-05
25 2.0991695E+00-8.1722E-02-8.1726E-02-0.0000E+00-8.4241E-02 3.0846E-01 2.6087E-05
26 2.0914340E+00-8.3724E-02-8.3729E-02-0.0000E+00-8.8341E-02 3.1329E-01 2.6610E-05
27 2.0838355E+00-8.5469E-02-8.5474E-02-0.0000E+00-9.2344E-02 3.1699E-01 2.6932E-05
28 2.0763925E+00-8.6937E-02-8.6941E-02-0.0000E+00-9.6263E-02 3.1946E-01 2.6998E-05
29 2.0691197E+00-8.8105E-02-8.8109E-02-0.0000E+00-9.9948E-02 3.2072E-01 2.6887E-05
30 2.0620282E+00-8.8951E-02-8.8956E-02-0.0000E+00-1.0339E-01 3.2070E-01 2.5076E-05
31 2.0551252E+00-8.9452E-02-8.9457E-02-0.0000E+00-1.0642E-01 3.1940E-01 2.5023E-05
32 2.0484146E+00-8.9583E-02-8.9588E-02-0.0000E+00-1.0900E-01 3.1676E-01 2.3536E-05
33 2.0418972E+00-8.9318E-02-8.9321E-02-0.0000E+00-1.1094E-01 3.1280E-01 2.1147E-05
34 2.0355709E+00-8.8626E-02-8.8630E-02-0.0000E+00-1.1216E-01 3.0748E-01 1.9440E-05
35 2.0294307E+00-8.7480E-02-8.7482E-02-0.0000E+00-1.1245E-01 3.0081E-01 1.5503E-05
36 2.0234697E+00-8.5844E-02-8.5846E-02-0.0000E+00-1.1171E-01 2.9276E-01 1.1283E-05
37 2.0176789E+00-8.3682E-02-8.3683E-02-0.0000E+00-1.0970E-01 2.8335E-01 6.6979E-06
38 2.0120481E+00-8.0953E-02-8.0954E-02-0.0000E+00-1.0631E-01 2.7251E-01 3.7512E-06
39 2.0065666E+00-7.7613E-02-7.7613E-02-0.0000E+00-1.0131E-01 2.6022E-01-4.0259E-06
40 2.0012241E+00-7.3607E-02-7.3606E-02-0.0000E+00-9.4618E-02 2.4636E-01-1.0741E-05
41 1.9960123E+00-6.8874E-02-6.8871E-02-0.0000E+00-8.6104E-02 2.3076E-01-1.8198E-05
42 1.9909271E+00-6.3336E-02-6.3333E-02-0.0000E+00-7.5800E-02 2.1314E-01-2.7663E-05
43 1.9859713E+00-5.6898E-02-5.6893E-02-0.0000E+00-6.3797E-02 1.9307E-01-4.2704E-05
44 1.9811595E+00-4.9431E-02-4.9424E-02-0.0000E+00-5.0426E-02 1.6984E-01-6.6669E-05
45 1.9765251E+00-4.0744E-02-4.0735E-02-0.0000E+00-3.6219E-02 1.4242E-01-1.0663E-04
46 1.9721292E+00-3.0521E-02-3.0509E-02-0.0000E+00-2.2083E-02 1.0907E-01-1.9366E-04
47 1.9680927E+00-1.8064E-02-1.8047E-02-0.0000E+00-9.3479E-03 6.6335E-02-4.4431E-04
48 1.9645663E+00 0.0000E+00 0.0000E+00-0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
The midplane positions of the plasma vacuum interface (PVI) are output. Here the vacuum construction can often fail because the conducting wall has not been sufficiently well built. The code will then stop and print and error message. For NOWALL=1 make sure the value of GWALL lies between the interior and exterior PVI values. If the pseudo-surfaces intersect on surfaces ~1-3, increasing the values of PVAC,QVAC can help. If the intersection occurs for values close to IVAC, decrease PVAC, QVAC.
PLASMA-VACUUM INTERFACE AND WALL EXTREMUM POSITIONS
INTERIOR WALL|INTERIOR PVI|EXTERIOR PVI|EXTERIOR WALL| WALL / PVI
4.02483E+00 4.16864E+00 6.15511E+00 6.27756E+00 1.13404E+00
***** BOPHYS, RECONSTRUCTED EQUILIBRIUM *****
I VP(I) PPI(I) PVPI(I) PARPVI(I) CJPI(I) CIPI(I) EQUI
1 2.2527E+00 4.1877E-03-2.1561E-03-0.0000E+00-7.8124E-03-2.0656E-02 8.8773E-01
2 2.2487E+00-4.3983E-03-4.9130E-03-0.0000E+00-1.2256E-02 1.6792E-02 5.5269E-02
3 2.2445E+00-7.7679E-03-7.9723E-03-0.0000E+00-1.5910E-02 3.0818E-02 1.2988E-02
4 2.2403E+00-1.1133E-02-1.1236E-02-0.0000E+00-1.9107E-02 4.4823E-02 4.6174E-03
5 2.2359E+00-1.4594E-02-1.4651E-02-0.0000E+00-2.1963E-02 5.9245E-02 1.9337E-03
6 2.2313E+00-1.8151E-02-1.8182E-02-0.0000E+00-2.4708E-02 7.4014E-02 8.3407E-04
7 2.2265E+00-2.1787E-02-2.1801E-02-0.0000E+00-2.7320E-02 8.9059E-02 3.2322E-04
8 2.2214E+00-2.5482E-02-2.5485E-02-0.0000E+00-2.9938E-02 1.0425E-01 6.6415E-05
9 2.2161E+00-2.9220E-02-2.9216E-02-0.0000E+00-3.2487E-02 1.1953E-01-6.1027E-05
10 2.2104E+00-3.2982E-02-3.2974E-02-0.0000E+00-3.5083E-02 1.3476E-01-1.1929E-04
11 2.2045E+00-3.6753E-02-3.6743E-02-0.0000E+00-3.7641E-02 1.4992E-01-1.4212E-04
12 2.1982E+00-4.0516E-02-4.0505E-02-0.0000E+00-4.0277E-02 1.6486E-01-1.3646E-04
13 2.1916E+00-4.4257E-02-4.4246E-02-0.0000E+00-4.2913E-02 1.7956E-01-1.1467E-04
14 2.1847E+00-4.7959E-02-4.7950E-02-0.0000E+00-4.5667E-02 1.9390E-01-9.1171E-05
15 2.1776E+00-5.1608E-02-5.1602E-02-0.0000E+00-4.8468E-02 2.0783E-01-6.3164E-05
16 2.1702E+00-5.5190E-02-5.5185E-02-0.0000E+00-5.1432E-02 2.2125E-01-4.1432E-05
17 2.1626E+00-5.8687E-02-5.8685E-02-0.0000E+00-5.4484E-02 2.3412E-01-2.2655E-05
18 2.1548E+00-6.2086E-02-6.2085E-02-0.0000E+00-5.7732E-02 2.4633E-01-5.9231E-06
19 2.1470E+00-6.5369E-02-6.5370E-02-0.0000E+00-6.1093E-02 2.5785E-01 5.9673E-06
20 2.1390E+00-6.8522E-02-6.8523E-02-0.0000E+00-6.4659E-02 2.6857E-01 1.2815E-05
21 2.1310E+00-7.1526E-02-7.1529E-02-0.0000E+00-6.8334E-02 2.7847E-01 1.9523E-05
22 2.1230E+00-7.4366E-02-7.4370E-02-0.0000E+00-7.2191E-02 2.8745E-01 2.2697E-05
23 2.1150E+00-7.7024E-02-7.7028E-02-0.0000E+00-7.6122E-02 2.9549E-01 2.4506E-05
24 2.1070E+00-7.9482E-02-7.9487E-02-0.0000E+00-8.0180E-02 3.0250E-01 2.6894E-05
25 2.0992E+00-8.1722E-02-8.1726E-02-0.0000E+00-8.4241E-02 3.0846E-01 2.6088E-05
26 2.0914E+00-8.3724E-02-8.3729E-02-0.0000E+00-8.8341E-02 3.1329E-01 2.6610E-05
27 2.0838E+00-8.5469E-02-8.5474E-02-0.0000E+00-9.2344E-02 3.1699E-01 2.6935E-05
28 2.0764E+00-8.6937E-02-8.6941E-02-0.0000E+00-9.6263E-02 3.1946E-01 2.6997E-05
29 2.0691E+00-8.8105E-02-8.8109E-02-0.0000E+00-9.9948E-02 3.2072E-01 2.6885E-05
30 2.0620E+00-8.8951E-02-8.8956E-02-0.0000E+00-1.0339E-01 3.2070E-01 2.5078E-05
31 2.0551E+00-8.9452E-02-8.9457E-02-0.0000E+00-1.0642E-01 3.1940E-01 2.5022E-05
32 2.0484E+00-8.9583E-02-8.9588E-02-0.0000E+00-1.0900E-01 3.1676E-01 2.3538E-05
33 2.0419E+00-8.9318E-02-8.9321E-02-0.0000E+00-1.1094E-01 3.1280E-01 2.1148E-05
34 2.0356E+00-8.8626E-02-8.8630E-02-0.0000E+00-1.1216E-01 3.0748E-01 1.9439E-05
35 2.0294E+00-8.7480E-02-8.7482E-02-0.0000E+00-1.1245E-01 3.0081E-01 1.5502E-05
36 2.0235E+00-8.5844E-02-8.5846E-02-0.0000E+00-1.1171E-01 2.9276E-01 1.1283E-05
37 2.0177E+00-8.3682E-02-8.3683E-02-0.0000E+00-1.0970E-01 2.8335E-01 6.6985E-06
38 2.0120E+00-8.0953E-02-8.0954E-02-0.0000E+00-1.0631E-01 2.7251E-01 3.7514E-06
39 2.0066E+00-7.7613E-02-7.7613E-02-0.0000E+00-1.0131E-01 2.6022E-01-4.0249E-06
40 2.0012E+00-7.3607E-02-7.3606E-02-0.0000E+00-9.4618E-02 2.4636E-01-1.0739E-05
41 1.9960E+00-6.8874E-02-6.8871E-02-0.0000E+00-8.6104E-02 2.3076E-01-1.8198E-05
42 1.9909E+00-6.3336E-02-6.3333E-02-0.0000E+00-7.5800E-02 2.1314E-01-2.7667E-05
43 1.9860E+00-5.6898E-02-5.6893E-02-0.0000E+00-6.3797E-02 1.9307E-01-4.2702E-05
44 1.9812E+00-4.9431E-02-4.9424E-02-0.0000E+00-5.0426E-02 1.6984E-01-6.6666E-05
45 1.9765E+00-4.0744E-02-4.0735E-02-0.0000E+00-3.6219E-02 1.4242E-01-1.0664E-04
46 1.9721E+00-3.0521E-02-3.0509E-02-0.0000E+00-2.2083E-02 1.0907E-01-1.9366E-04
47 1.9681E+00-1.8064E-02-1.8047E-02-0.0000E+00-9.3479E-03 6.6335E-02-4.4431E-04
48 1.9646E+00 0.0000E+00 0.0000E+00-0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
PLASMA VOLUME = 8.316995E+01 m^3
AFVEN FREQUENCY SQUARED x MU_0 x MASS DENSITY = 2.827659E+01
1
####### TERPSICHORE ######## STABIN ####### TERPSICHORE ######
NSTA MMS NSMIN NSMAX MODEL
3 25 -20 22 0
TABLE OF R AND Z COEFFICIENTS FOR STABILITY
M= 0 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 ...N
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1-20
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-19
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-18
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1-17
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-16
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-15
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1-14
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-13
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-12
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1-11
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-10
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -9
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -8
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -7
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -5
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -4
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -3
0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -2
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 4
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 6
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 7
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 8
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 9
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 10
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 11
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 13
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 14
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 15
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 16
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 17
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 18
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 19
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 21
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 22
NEV NITMAX AL0 EPSPAM IGREEN MPINIT
1 4500 -5.000E-05 1.000E-04 0 0
NON-INTERACTING ENERGETIC PARTICLE MODEL WITH REDUCED KINETIC ENERGY: modelk = 0
Number of Modes in the Stability Table: 3821463
* FILTER
MATRICES FROM TERPSICHORE ARE USED FOR PAMERA
A -> -A ETC BECAUSE OF SIGN(BJAC)
* FILTER
382 60 60 61
382 60 382
1
EIGENVALUE SHIFT = -5.00000E-05
RAYLEIGH QUOTIENT = 4.59014E-07
EIGENVALUE FROM NORMALIZATION = 4.58935E-07
MAX. NUMBER OF ITERATIONS = 4500
ITERATIONS DONE = 69
NUMBER OF NEGATIVE EIGENVALUES = 0
NUMBER OF NEGATIVE X - PIVOTS = 0
NUMBER OF NEGATIVE Y - PIVOTS = 0
NUMBER OF NON CONVERGED X COMPONENTS = 0
NUMBER OF NON CONVERGED Y COMPONENTS = 0
CONVERGED IF ERROR SMALLER THEN 1.00000E-04
1
WP = 8.79030E-07 WK = 1.91504E+00 EIGENVALUE FROM WP / WK = 4.59013931E-07
GROWTH RATE = -1.27408797E-04
I Q(SI) DELTAW(SI) MERCIER CRITERION
KO-standard KO-lambda NON-INTER_st NON-INTER_lamba
1 4.43812829E+00 2.37746375E-10 -6.8703E+00 -9.5868E-01 -6.8695E+00 -9.5868E-01
2 4.22000941E+00 2.37947000E-10 2.7632E-01 2.3363E+02 2.7632E-01 2.3364E+02
3 4.03693189E+00 2.12455365E-10 1.0749E+00 3.3922E+00 1.0749E+00 3.3923E+00
4 3.88289834E+00 1.95350781E-10 1.2067E+00 2.1034E+00 1.2068E+00 2.1034E+00
5 3.74995485E+00 1.83122841E-10 1.2945E+00 1.4915E+00 1.2946E+00 1.4916E+00
6 3.63353320E+00 1.73508441E-10 1.3606E+00 1.1358E+00 1.3607E+00 1.1359E+00
7 3.52948521E+00 1.65126083E-10 1.4172E+00 9.1233E-01 1.4173E+00 9.1245E-01
8 3.43565310E+00 1.57211457E-10 1.4785E+00 7.6657E-01 1.4788E+00 7.6671E-01
9 3.34995762E+00 1.49630907E-10 1.5549E+00 6.7094E-01 1.5553E+00 6.7109E-01
10 3.27138305E+00 1.42661504E-10 1.6548E+00 6.0869E-01 1.6553E+00 6.0887E-01
11 3.19876580E+00 1.37077176E-10 1.7816E+00 5.6981E-01 1.7822E+00 5.7001E-01
12 3.13152492E+00 1.33945692E-10 1.9357E+00 5.4717E-01 1.9365E+00 5.4740E-01
13 3.06885404E+00 1.34773375E-10 2.1126E+00 5.3622E-01 2.1135E+00 5.3648E-01
14 3.01029867E+00 1.44468610E-10 2.3035E+00 5.3308E-01 2.3047E+00 5.3336E-01
15 2.95518607E+00 1.60206664E-10 2.4966E+00 5.3575E-01 2.4980E+00 5.3607E-01
16 2.90309176E+00 1.72374210E-10 2.6787E+00 5.4273E-01 2.6804E+00 5.4307E-01
17 2.85341890E+00 1.82118477E-10 2.8336E+00 5.5263E-01 2.8355E+00 5.5301E-01
18 2.80577732E+00 1.90278222E-10 2.9478E+00 5.6447E-01 2.9499E+00 5.6489E-01
19 2.75965999E+00 1.95838771E-10 3.0119E+00 5.7778E-01 3.0142E+00 5.7824E-01
20 2.71475010E+00 1.98693740E-10 3.0224E+00 5.9190E-01 3.0249E+00 5.9239E-01
21 2.67066248E+00 2.01136780E-10 2.9795E+00 6.0579E-01 2.9822E+00 6.0633E-01
22 2.62719029E+00 2.02997209E-10 2.8927E+00 6.1963E-01 2.8954E+00 6.2021E-01
23 2.58409194E+00 2.03780862E-10 2.7749E+00 6.3441E-01 2.7776E+00 6.3504E-01
24 2.54128717E+00 2.02960663E-10 2.6367E+00 6.4955E-01 2.6394E+00 6.5022E-01
25 2.49867860E+00 1.98965926E-10 2.4881E+00 6.6481E-01 2.4908E+00 6.6554E-01
26 2.45630701E+00 1.98401563E-10 2.3382E+00 6.8003E-01 2.3408E+00 6.8082E-01
27 2.41420158E+00 2.05213416E-10 2.1938E+00 6.9528E-01 2.1964E+00 6.9612E-01
28 2.37250300E+00 2.17281136E-10 2.0601E+00 7.1055E-01 2.0627E+00 7.1145E-01
29 2.33133873E+00 2.36233313E-10 1.9405E+00 7.2600E-01 1.9430E+00 7.2696E-01
30 2.29092164E+00 2.56981703E-10 1.8383E+00 7.4193E-01 1.8407E+00 7.4295E-01
31 2.25144768E+00 2.94941689E-10 1.7551E+00 7.5854E-01 1.7575E+00 7.5962E-01
32 2.21317612E+00 3.55930337E-10 1.6922E+00 7.7593E-01 1.6947E+00 7.7707E-01
33 2.17634777E+00 4.06176425E-10 1.6519E+00 7.9456E-01 1.6543E+00 7.9577E-01
34 2.14125117E+00 4.48949280E-10 1.6370E+00 8.1459E-01 1.6395E+00 8.1587E-01
35 2.10815869E+00 4.94140712E-10 1.6523E+00 8.3554E-01 1.6549E+00 8.3688E-01
36 2.07738124E+00 5.22929601E-10 1.7059E+00 8.5607E-01 1.7086E+00 8.5748E-01
37 2.04921858E+00 4.92256140E-10 1.8140E+00 8.7409E-01 1.8170E+00 8.7557E-01
38 2.02400351E+00 2.59462427E-10 2.0082E+00 8.8566E-01 2.0116E+00 8.8719E-01
39 2.00206176E+00 -9.72057817E-10 2.3391E+00 8.7912E-01 2.3431E+00 8.8068E-01
40 1.98374377E+00 -4.67168332E-09 2.9609E+00 8.6188E-01 2.9662E+00 8.6346E-01
41 1.96938776E+00 -8.93641353E-09 4.3885E+00 8.5530E-01 4.3966E+00 8.5690E-01
42 1.95933334E+00 -1.30827525E-08 8.9628E+00 8.6807E-01 8.9796E+00 8.6971E-01
43 1.95387564E+00 -1.63704400E-08 4.7993E+01 9.0686E-01 4.8081E+01 9.0853E-01
44 1.95324596E+00 -1.70716143E-08 1.1041E+02 9.8769E-01 1.1060E+02 9.8940E-01
45 1.95752895E+00 -1.00222728E-08 6.9384E+00 1.1485E+00 6.9488E+00 1.1503E+00
46 1.96658102E+00 2.74350413E-08 2.0854E+00 1.4798E+00 2.0877E+00 1.4815E+00
47 1.97980793E+00 3.62984082E-07 9.5738E-01 2.3119E+00 9.5801E-01 2.3136E+00
48 1.99617566E+00 3.93353419E-07 4.7862E-01 5.3743E+00 4.7876E-01 5.3767E+00
49 2.00435953E+00 5.81525833E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
50 2.00435953E+00 5.64196009E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
51 2.00435953E+00 5.32722605E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
52 2.00435953E+00 4.91835452E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
53 2.00435953E+00 4.46713465E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
54 2.00435953E+00 4.01811090E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
55 2.00435953E+00 3.60269666E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
56 2.00435953E+00 3.23906982E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
57 2.00435953E+00 2.93525053E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
58 2.00435953E+00 2.69309125E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
59 2.00435953E+00 2.51186683E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
60 2.00435953E+00 2.39114176E-09 0.0000E+00 0.0000E+00 0.0000E+00 0.0000E+00
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
125 1 -5 5.753E-07 5.620E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
126 2 -5 6.936E-07 1.146E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
127 3 -5 1.533E-06 5.141E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
128 4 -5 8.595E-07 1.806E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
129 5 -5 7.744E-07 1.232E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
150 1 -2 4.756E-06 5.622E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
151 2 -2 1.128E-05 3.652E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
152 3 -2 4.935E-06 1.086E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
153 4 -2 8.457E-06 1.455E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
154 5 -2 2.328E-06 6.603E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
155 6 -2 4.414E-06 2.303E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
156 7 -2 2.667E-06 2.842E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
157 8 -2 9.430E-07 1.522E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
158 9 -2 1.384E-06 1.525E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
159 10 -2 1.034E-06 6.997E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
175 0 1 3.048E-07 9.203E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
176 1 1 6.173E-05 4.819E-04
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
177 2 1 5.712E-02 9.656E-01
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
178 3 1 2.386E-04 4.788E-04
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
179 4 1 4.177E-05 1.381E-04
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
180 5 1 2.092E-05 7.721E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
181 6 1 2.621E-05 3.647E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
182 7 1 2.968E-05 1.855E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
183 8 1 9.848E-06 6.348E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
184 9 1 7.264E-06 3.249E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
185 10 1 5.858E-06 5.833E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
186 11 1 3.577E-06 4.361E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
187 12 1 3.908E-06 1.907E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
188 13 1 2.612E-06 1.533E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
189 14 1 1.278E-06 7.276E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
190 15 1 1.185E-06 5.040E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
201 0 4 1.389E-07 1.379E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
202 1 4 1.761E-06 2.961E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
203 2 4 6.392E-06 1.574E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
204 3 4 1.139E-05 4.007E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
205 4 4 3.429E-05 7.867E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
206 5 4 2.263E-05 5.204E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
207 6 4 1.659E-05 1.032E-04
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
208 7 4 1.767E-04 1.409E-04
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
209 8 4 1.061E-02 5.033E-02
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
210 9 4 4.693E-04 1.192E-04
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
211 10 4 1.321E-04 4.996E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
212 11 4 4.378E-05 3.913E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
213 12 4 3.860E-05 1.682E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
214 13 4 3.567E-05 9.977E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
215 14 4 1.431E-05 4.416E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
216 15 4 6.230E-06 1.762E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
217 16 4 3.597E-06 2.953E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
218 17 4 1.939E-06 1.453E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
219 18 4 1.932E-06 6.095E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
229 2 7 4.978E-07 2.815E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
230 3 7 1.020E-06 1.986E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
231 4 7 2.645E-06 5.058E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
232 5 7 4.877E-06 1.025E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
233 6 7 1.267E-05 1.904E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
234 7 7 1.244E-05 1.464E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
235 8 7 6.037E-06 3.859E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
236 9 7 3.449E-05 1.733E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
237 10 7 2.670E-05 2.056E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
238 11 7 3.279E-05 1.974E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
239 12 7 6.098E-05 3.725E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
240 13 7 2.453E-04 1.317E-04
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
241 14 7 5.414E-03 1.502E-02
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
242 15 7 4.143E-04 5.581E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
243 16 7 1.820E-04 4.842E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
244 17 7 8.016E-05 2.648E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
245 18 7 2.853E-05 1.579E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
246 19 7 2.071E-05 3.511E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
247 20 7 6.627E-06 1.213E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
248 21 7 2.784E-06 1.168E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
249 22 7 1.073E-06 3.921E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
259 6 10 2.188E-06 3.035E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
260 7 10 3.648E-06 4.861E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
261 8 10 5.329E-06 5.743E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
262 9 10 2.992E-06 3.474E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
263 10 10 2.216E-06 1.166E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
264 11 10 6.188E-06 4.794E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
265 12 10 6.382E-06 8.085E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
266 13 10 3.990E-06 4.016E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
267 14 10 4.813E-06 1.579E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
268 15 10 3.281E-05 9.210E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
269 16 10 2.280E-05 9.131E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
270 17 10 2.787E-05 1.124E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
271 18 10 5.326E-05 3.733E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
272 19 10 1.817E-04 7.892E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
273 20 10 2.060E-03 4.325E-03
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
274 21 10 2.444E-04 2.840E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
275 22 10 1.404E-04 2.770E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
276 23 10 6.422E-05 1.505E-05
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
277 24 10 2.033E-05 6.874E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
278 25 10 5.391E-06 2.303E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
287 8 13 8.347E-07 1.054E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
288 9 13 1.610E-06 1.728E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
289 10 13 1.882E-06 2.071E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
290 11 13 6.034E-07 1.626E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
291 12 13 1.408E-06 5.133E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
292 13 13 2.184E-06 2.607E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
293 14 13 2.477E-06 2.844E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
294 15 13 8.666E-07 1.577E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
295 16 13 2.464E-06 6.471E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
296 17 13 3.427E-06 3.168E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
297 18 13 1.670E-06 3.034E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
298 19 13 7.473E-07 1.511E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
299 20 13 8.261E-06 2.818E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
300 21 13 2.150E-05 4.045E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
301 22 13 1.036E-05 4.978E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
302 23 13 4.491E-06 5.259E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
303 24 13 4.098E-06 5.872E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
304 25 13 1.211E-05 7.581E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
321 16 16 7.267E-07 1.363E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
322 17 16 1.456E-06 6.030E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
323 18 16 1.242E-06 4.573E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
324 19 16 7.798E-07 1.157E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
325 20 16 7.229E-07 1.262E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
326 21 16 1.820E-06 9.122E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
327 22 16 2.655E-06 5.256E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
328 23 16 1.763E-06 7.678E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
329 24 16 1.284E-06 1.247E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
330 25 16 2.774E-06 1.527E-06
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
354 23 19 1.344E-06 3.825E-07
PARVEC: LC, MS(LC), NS(LC), XIMAX, ETAMAX
355 24 19 1.094E-06 1.865E-07
Modes in the Stability Table in which Xi,Eta Fourier Amplitude exceeds:
1.00000E-06 121
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-20
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-19
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-18
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-17
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-16
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-15
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-14
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-13
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-12
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-11
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0-10
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -9
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -8
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -7
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -6
0 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -5
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -4
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -3
0 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -2
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 6
0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 7
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 8
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 9
0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 10
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 11
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12
0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 13
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 14
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 15
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 16
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 17
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 18
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 19
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 21
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 22
Timing in subroutine eqinvm = 7.199300E-02 seconds
Timing in subroutine veqrec = 7.601691E+01 seconds
Timing in subroutine mtaskb = 1.216520E+02 seconds
Timing in subroutine mtaska = 1.029535E+02 seconds
Timing in solver subroutines= 5.071488E+01 seconds
Total running time = 3.514094E+02 seconds
2 1 1.68816069E-02
8 4 4.41167945E-03
14 7 2.46410908E-03
20 10 1.15232245E-03
9 4 4.51602990E-04